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探测 器 VLBI 信号 仿真 方法 : 
RBZ, WAR, $51, EWI, KAI, WA! 
(1 中 国 科学 院 上 海天 文 台 ， 上 海 200030; 2, 中国 科学 院 射电 天 文 重点 实验 室 ， 江 苏 “ 南京 
210008; 3. 上 海 市 导航 定位 重点 实验 室 ， 上 海 200030; 4. 中 国 科学 院 大 学 ， 北 京 100049) 


摘要 : AAR VLBI 对 探测 器 的 数据 处 理 方 法 并 评估 其 处 理 能 力 ， 现 有 做 法 需 基 于 测 站 
VLBI 终端 接收 的 探测 器 信号 进行 。 采 用 信号 仿真 方法 可 以 根据 设计 轨道 与 信 标 特点 利用 计 
算 机 生成 需要 的 探测 器 VLB1 信号 ， 相 比试 验 观测 具有 独特 优势 。 观 测试 验 时 ， ees 
收 探 测 器 下 行 射频 信号 ， 通 过 混 频 变换 为 中 频 信 号 ， 再 进行 数字 化 采集 数据 。 仿 真 信号 

为 减少 计算 规模 ,直接 构造 数字 中 频 信号 ,二 次 采样 后 提取 通道 信号 eras 
仿真 信号 。 由 于 探测 器 相对 VLBI 测 站 运动 ， 测 站 接收 的 探测 器 下 行 信号 反映 目标 的 视 向 速 

度 变 化 。 根 据 探测 器 信号 时 延 模型 ， 研 究 了 探测 器 点 频 信 号 和 有 限 带宽 信号 多 普 勒 效应 的 模 
拟 方 法 。 通 过 仿真 信号 与 VLB1 终端 接收 信号 的 对 比 以 及 对 VLBI 相关 处 理 结果 的 分 析 ， 验 证 
了 探测 器 VLBI 信号 仿真 方法 的 可 行 性 , 为 后 续 仿真 技术 在 月 球 与 火星 探测 器 中 的 VLBI1 测定 
轨 应 用 黄 定 了 技术 基础 。 


关键 字 : 信号 仿真 ; 深 空 探 测 ; VLBI; 信 标 
中 图 分 类 号 : 文献 标识 码 : A 文章 编号 : 1672-7673 


其 长 基线 干涉 测量 由 (Very Long Baseline Interferometry, VLBI) 是 20 世纪 60 年 代 后 期 
出 现 的 一 种 射电 天 文 技术 。 它 具有 空间 分 辨 率 高 、 测 量 精 度 高 的 特点 ， 现 已 广泛 应 用 于 天 体 
物理 、 天 体 测 量 、 天 文 地 球 动力 学 和 深 空 探测 领域 名。 我 国 于 2007 年 首次 在 探 月 工程 中 使 
H VLBI R, Æ T HE VLBI A (Chinese VLBI Network, CVN) 包括 上 海 余 山 站 、 北 京 密 
云 站 、 云 南 昆明 站 、 乌 鲁 木 齐 南山 站 射电 望远镜 和 上 海 VLBI 数据 处 理 中 心 。 中 国 VLBI 网 已 
参加 了 控 月 工程 嫦娥 系列 探测 器 的 多 次 任务 , 对 关键 弧 段 的 快速 精密 测定 轨 发 挥 了 不 可 替代 
的 重要 作用 铝 。 

VLBI 测 站 由 天 线 、 接 收 机 、 数 据 采 集 终端 、 时 间 频 率 系 统 等 组 成 ， 在 工作 时 ， 各 测 站 接 
收 机 天 线 汇集 的 射频 信和 号 放大 后 ,经 独立 本 振 进 行 混 频 得 到 中 频 信号 ， 然 后 传送 至 数字 化 终 
端 ， 高 速 模 数 采 样 后 作 数 字 基带 转换 ， 得 到 所 需 频 率 通道 的 基 频 信号 ， 最 后 经 过 编码 ， 记 录 
在 磁盘 或 者 直接 通过 网 络 传送 到 VLBI 中 心口。 

为 研究 VLBI 中 心 对 不 同 探测 器 的 测量 方法 ， 评 估 分 析 中 心 各 功能 模块 的 能 力 及 整个 系 
统 的 正确 性 、 可 靠 性 ,目前 的 常规 做 法 是 在 每 次 飞行 任务 前 开展 大 量 实际 目标 联 测 。 对 单 目 
标 探测 器 的 实际 观测 演练 可 评估 测 站 与 系统 的 可 靠 性 。 但 对 于 那些 全 新 的 深 空 探测 任务 , 在 
ene ieee era nid lances: 比如 在 后 续 探 月 工程 中 , 将 要 开展 针 
对 月 球 轨道 的 交会 对 接 ， 届 时 VLBI 需要 进行 对 月 球 距离 的 两 个 动态 目标 的 同 波 束 测量 。 显 
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号 仿真 方法 ， 则 可 根据 探测 器 的 数量 、 轨 道 


、 信 标 等 特点 ， 模 拟 生成 VLBI 终端 接收 的 不 同 


探测 器 信号 ， 提 供 月 球 轨道 双 目 标 同 波束 VLBI 测量 数据 ， 用 于 研究 和 改进 VLBI 处 理 


模块 。 


通过 


寺 仿 真 与 实际 


演练 两 种 方法 结合 ， 可 较 好 地 记 


EF 估 VLBI 中 心 与 整个 系统 对 双 CE) 目标 探 


测 器 测量 能 力 。 因 此 ， 研 究 VLBI 信和 号 仿真 技术 有 重要 意义 与 迫切 工程 需求 。 
在 探 月 工程 中 ， 探 测 器 下 行 VLBI 射频 信号 为 S/X 频段 ， 采 用 调 相 方式 ， 调 制 信号 在 形 
式 上 主要 包括 点 频 信 号 和 有 限 带宽 信号 两 类 。 在 仿真 信号 时 ， 没 有 直接 仿真 射频 信号 ， 而 是 


先 将 信号 从 射频 变换 至 数字 中 频 ， 然 后 再 对 中 频 信号 采样 ， 提 取 通 道 信 号 ， 最 终生 成 用 于 验 
证 VLBI 数据 处 理 系统 的 仿真 数据 ， 这 样 既 完整 地 保留 了 调制 信号 频谱 主要 特性 ， 又 避免 了 


直接 处 理 需 极 高 采样 率 的 射频 信号 。 
探测 器 信号 结构 
在 探 月 工程 中 ， 探 测 器 VLB 下 行 射频 信号 采用 调 相 方式 ， 包 括 遥 测 信号 、 单 向 差分 测 
BE (DOR) 信号 、 测 距 信号 、 残 留 遥 控 信 号 〈 图 1)。 其 中 ， 下 行 载波 、 测 距 信号 〈 测 距 主 音 
和 测 距 次 音 ) 和 单 向 差分 测 距 信 号 为 点 频 信 号 ; 遥测 信号 和 遥控 信号 为 经 过 相位 调制 的 有 限 
带宽 信号 。 这 些 信号 统一 调制 在 下 行 信号 的 主 载波 上 ， 数 学 模型 外 为 
s(t) = J2P sin@aft + Od(t) + 0,d,(t) 


+ gsin(2zft) + hsin(2nt,t)) 
其 中 ，:Pr 为 总 功率 ; f, 为 载波 频率 ，di(7) 为 遥测 信号 ，0 为 遥测 信号 的 调制 指数 ，c(z) 
为 遥控 信号 ; @ 为 膛 控 信号 的 调制 指数 ， 为 测 距 信号 ;9 为 测 距 信号 的 调制 指数 ， 石 为 


DOR 信号 ; o, 为 DOR 信号 的 调制 指数 。 
下 行 
载波 | AA WE ,， 测 距 测 距 


DOR 音 1 DOR 音 2 


遥控 KÄ 


8kHz 16~ 65.536 100kHz 500kHz 3.8MHz 19.2MHz  f-fc 
20kHz kHz 
Al xX 波段 下 行 信号 上 边 带 频谱 示意 图 
Fig. The spectrum of X band downlink signal 


2 仿真 原理 与 实现 方法 


2.1 探测 器 信号 时 延 模 型 构造 
VLBI nisi sea > LATERM, 因此 在 仿真 时 要 考虑 相对 运动 效应 引起 的 测 站 接 
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收 探测 器 信号 的 时 间 延 退 〈 时 延 ) 及 其 变化 加。 时 延 包 含 几何 时 延 、 引 力 时 延 、 天 线 时 延 修 
正 、 源 结构 时 延 修正 、 中 性 大 气 与 电离 层 时 延 修正 等 四， 其 中 ， 几 何 时 延 是 理论 时 延 的 主要 


部 分 ， 随 目标 的 运动 而 变化 。 


设 地 心 0、 测 站 4 及 探测 器 8B (图 2)。g0(t) 表 示 t 时 刻 参 考点 (地 心 ) 接收 到 信号 波 


前 的 光 行 时 ， gO 表示 tt 时 刻 参 考点 〈 地 心 ) 信号 波 前 超前 到 达 测 站 的 光 行 时 ， g 0) K 


示 志 时 刻 参考 点 〈 测 站 ) 接收 到 信号 波 前 的 光 行 时 。 若 时 刻 从 探测 器 8 发 出 信号 ， 同 一 
波 前 到 达 地 心 0、 测 站 A NAIA AC, 如 ， 则 根据 定义 eo (ty). 2: (ty) Raey(t,) BAA 


golto) = to — by 


gto) =t,-t 
g(t) = -tl 
显然 ， 


(2) 


E(t) = golto) + gi(to) (3) 


图 2 VLBI 时 延 示意 图 


Fig. 2 The diagram of VLBI 


g(t) = b + bt + bt? + Bt? 
KH, by, b. b AY REALS ARB. 


delay 


在 实际 应 用 中 ， 每 个 观测 弧 段 的 时 延 模 型 都 以 五 次 时 延 多 项 式 表示 : 


+ b,t* + bt, (4) 


根据 实际 时 延 模型 g0(t), g1(z) ,用 最 小 二 乘法 作 五 次 多 项 式 拟 合 , 得 时 延 模 型 g) 


2.2 时 延 模型 线性 近似 


为 减少 计算 量 ， 在 短 时 间 内 假设 时 延 模型 随时 间 线 性 变化 ， 即 


gt) = bh + bt 


， (5) 


Ki, ABD, bi 在 短 时 间 内 为 恒 常 数 ， 忽 略 时 延 模 型 8,(t) 高 次 项 。 


在 短 时 间 内 用 (5)〉 式 表示 线性 化 时 延 模 型 存在 


的 时 间 长 度 和 五 次 时 延 多 项 式 中 高 次 项 系数 。 相 关 处 至 


3 


定 误差 ， 误差 大 小 取决 于 用 线性 近似 
E 机 "使 用 皮 秒 量 级 的 时 延 模型 精度 足 


HAFIJ 


hinaXiv 合 作 其 


90 ”以 满足 应 用 需求 , 故 误差 值 在 0.01ps (1e-14s) 以 内 的 模型 误差 即 可 视 为 满足 仿真 精度 要 求 。 
91 ”在 较 短 的 观测 时 间 内 , 时 延 模型 g,(t) 的 二 次 项 是 线性 近似 的 主要 误差 项 , HAMS, 一 般 为 
92 ”e-11 量 级 , 要 想 使 得 该 项 引起 的 时 延误 差 小 于 0.01ps, 选 定 的 线性 化 时 间 长 度 应 小 于 0.01s。 
93 ”其 他 高 次 项 的 系数 乌 , b bs IRE e-14 量 级 ， 由 此 引起 的 误差 可 忽略 。 


94 
95 ”2.3 构造 中 频 信和 号 
96 探测 器 下 行 载波 信号 人 点 频 信 号 和 有 限 带 宽 信 号 。 其 中 点 频 信 号 包括 : 载 


97 ” 波 、 测 距 信 号 和 单 向 差分 测 距 信 号 。 有 限 带宽 信号 包括 遥测 信号 和 遥控 信号 ， 是 以 正弦 波 为 
98 ”载波 的 调 相 信号 ， 调 制 方 oT ee ane rit 
99 ”变化 传递 数字 信息 。 对 于 二 进 制 相 移 键 控 信号 ，0 以 0/ 元 交替 变化 表示 二 进 制 位 0/1， 1 
100 ”对 应 相位 0，-1 对 应 相位 元 。 因 此 遥测 和 遥控 信号 又 可 表示 为 4sin(2zxzz) ，4 随 相位 如 的 
101 ”变化 而 变化 。 

102 为 完整 地 保留 调制 信号 及 其 边 带 频谱 ， 选 择 中 频 信号 约 为 60MHz。 有 限 带 宽 信 号 的 码 
103 ”速率 最 大 为 2048bit/s， 码 元 变化 周期 是 采样 周期 的 两 万 九 千 多 倍 。 在 一 个 码 元 周期 内 ， 有 
104 ，” 限 带宽 信号 保持 不 变 ， 对 其 采样 可 仿真 遥测 和 遥控 有 限 带宽 信和 号。 因此， 需 选择 合适 的 时 间 


105 ”尺度 对 点 频 信 号 和 有 限 带 宽 信 号 采样 ， 仿 真 中 频 信号。 
106 为 简约 起 见 , 根据 (1) 式 探测 器 转发 的 下 行 信号 一 般 形式 ,分 析 其 中 典型 的 点 频 信和 号 、 
107 ”有 限 带 宽 信号 多 普 勒 效应 影响 。 信 号 表达 式 如 下 : 
~ 108 f(t) = sinzft + 0Asin aft) + hsin(2xt;t)), (6) 
= 109 Sh, £, 分 别 为 载波 和 点 频 单 向 差分 测 距 信 号 ，£ 为 有 限 带宽 遥测 信号 ，4 为 遥测 
I ono 。 信号 的 幅 值 ，Q, 办 为 调制 系数 。 
> 111 测 站 志 时 刻 接收 的 信号 A(z) 是 探测 器 在 上 一 g(t) 时 刻 发 出 的 ， 短 时 间 内 假设 时 延 模 
112 ”型 随时 间 线 性 变化 ， 用 g(t) Rg, (0), FO 与 探测 器 转发 的 了 (t) 关系 是 
三 113 F(t) = f(t - g,¢)), (7) 
© 114 将 (5) 式 和 (6) 式 代 入 (7) 式 ， 得 到 : 
115 F(t) = sinQxf((l — BE — h) + GAsin ZAU- Ak -= h). a 
+ h sin 2nf,((1 — b)t - b )) 
116 对 比 (6) 式 ， 测 站 接收 信号 EO 的 载波 和 信号 E, 已 发 生 频率 偏 移 ， 体 现 了 信号 传 
117 。 输 中 由 于 相对 运动 产生 的 多 普 勒 效应 。 ME ONIN LEME + E) b(f + f) 
118 根据 〈《8) 式 ， 探 测 器 下 行 信号 在 频 域 上 表现 为 关于 载波 左右 对 称 ， 其 频谱 主要 分 布 在 
119 [0 -4)(4-f£f,,),0-4)44+f,,)], HPL, Nmax(h, 6), b << 1. 
120 ”构造 中 频 信 号 了 (t) 如 下 : 
121 f(t) = nO) a 
+ h sin 2nf,((1 — &)t - b)) 
122 Sth, £ 为 信号 起 始 频 点 ( — b) (E - £,.)- 
123 对 比 C9) 式 与 (8) 式 ， 只 是 载波 频率 发 生变 化 ,其 余 信号 不 变 。 即 通过 降低 载波 频率 ， 
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将 信号 从 射频 搬移 到 中 频 ， 且 信号 A, 无 前 后 保持 不 变 。 
综 上 , VLBI 测 站 接收 信号 是 添加 时 延 的 射频 信号 , 由 (1) 式 添 加 时 延 并 降低 载波 频率 ， 
得 到 包含 多 普 勒 效应 的 中 频 信号 ， 形 式 如 下 : 
S(t) = J2P, sin(a(Al - 4) - FW -— £A) 
+ Asin 2zf,(0 — b)t - h) + 6A sin 2nf,( — b)t — b) (10) 
+ ¢, sin 22f((1 - bt - h) + ġ sin 2x6((1 -— bt - h )) 
th, A WEES, A ERA SRE. 


2.4 VLBI 终端 信号 仿真 

图 3 显示 了 仿真 VLBI 终端 接收 信号 实现 流程 。 先 仿真 带 有 多 普 勒 效应 的 探测 器 VLBI 下 
行 中 频 信号 ， 再 滤波 得 到 某 个 通道 的 频率 区 间 ， 通 过 和 希 尔 伯 特 变换 保留 正 频率 部 分 ， 再 通过 
数字 变频 移 频 至 基带 , 经 希 尔 伯 特 逆 变 换取 实 信号 , 最 后 通过 重 采 样 降低 采样 率 得 到 通道 实 
信和 号。 通道 宽度 为 2MHz 时 ， 采 样 率 取 4MHz。 


( 开始 ) 


生成 中 频 信号 


图 3 信号 仿真 流程 图 


Fig. 3 Steps of signal simulation 


3. 数值 试验 


3.1 仿真 探测 器 下 行 中 频 信号 

根据 上 述 方法 采用 某 月 球 探测 器 信 标 实际 参数 仿真 了 探测 器 VLBI 下 行 中 频 信 号 ， 其 中 
包含 测 距 信号 (500kHz 主音 频率 )、 遥 测 信号 《〈 副 载波 65.536kHz 且 码 速率 256bit/s) 和 两 
对 单 向 差分 测 距 信 号 。 高 频 信号 起 始 频 点 为 8460MHz， 降 低 载 波 频率 后 生成 探测 器 下 行 
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速 傅 里 时 变换 点 165000。 
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Fig. 4 The spectrum of 
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downlink probe quasi-intermediate-frequency signal 


图 4 显示 带宽 58.5MHz 的 探测 器 下 行 中 频 信号 频谱 。 图 中 ， 箭 头 3 为 载波 信号 ，2 和 4 
a 第 1 对 单 向 差分 测 距 信号 CDOR1); 1 和 5 是 距 载波 约 19.27MHz 的 
第 2 对 DOR 信号 (DOR2)， 载 波 两 侧 500kHz 处 是 测 距 主音 人 信号。 另外， 有限 带宽 的 遥测 信 
E 
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以 下 。 


3.2 PJ ÍT 


通过 仿真 信号 
号 仿真 方法 的 可 行 性 。 


测 器 VLBI 信 


与 VLBI 终端 接收 信号 


图 中 载波 两 侧 对 应 频 点 幅 值 产生 细微 差别 。 其 他 信号 谐 波 项 能 量 较 小 


>» i 


晶 度 基本 在 -180dB 


的 对 比 以 及 对 VLBI 相关 处 理 结果 的 分 析 ， 验 证 了 探 


图 5 中 (a) 为 实际 信号 单 通 道 自 功率 谱 ， 从 中 可 观察 到 0.96MHz 处 的 载波 信号 ; 距离 
载波 500kHz、65.536kHz 的 测 距 主 音信 号 和 遥测 信号 ; 其 他 为 信号 谐 波 项 。 
图 5 P (b) 为 仿真 VLBI 终端 接收 信号 自 功率 谱 (采样 率 4M, FFT 点 262144, 采样 1s )， 


包含 载波 信号 、 测 距 


1024bit/s)、 高 斯 白 


且 关 于 测 距 主 
仿真 信号 


音 对 称 的 为 遥 


4S (500kHz = 


与 对 应 的 实际 信 
配 与 实际 信号 一 致 ， 噪 声 信 和 号 的 功率 和 


Ul es 号 谐 波 项 。 


号 对 比 表 明 ， 


BA 波 信和 号 、 测 距 信 号 、 遥测 ao 
分布 也 与 实际 信号 相同 ， 验 证 了 仿真 方法 的 可 行 性 。 


音频 率 )、 遥 测 信号 〈 副 载波 65.536kHz 且 码 速率 
噪声 信号 。 频 率 0.968MHz 处 为 载波 信号 ， 测 距 主音 信 
测 信号 分 别 距 离 载 波 500kHz、65.536kHz, 关于 载波 对 称 。2 


号 和 有 限 带 宽 的 遥 


E 离 测 距 主音 65.536kHz、131.072kHz 


的 频率 和 功率 分 
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(a) (b) 
图 5 单 通道 自 功 率 谱 图 (A) 实际 信号 ; (B) 仿真 信号 


Fig. 5 The auto-power spectrum of channel signal. (A) The auto-power spectrum of real signal; 


(B) The auto-power spectrum of simulation signal 


探测 器 VLBI 下 行 中 频 信 号 分 别 添加 不 同时 延 模 型 ， 仿 真 两 台 站 终端 接收 信号 《采样 率 
4M, FFT 点 512， 采 样 1 秒 )。 经 过 VLBI 相关 处 理 ， 可 得 互 功率 谱 幅 度 和 相位 图 (图 6)。 图 
6 中 (C) 的 幅度 谱 显 示 ， 载 波 信号 和 遥测 信号 在 [0.8MHz,1.2MHz] 频 段 ， 遥 测 信号 的 谐 波 项 
和 测 距 信号 在 [0.4MHz,0.6MHz]，[1.4MHz,1.6MHz] 频 段 。 图 6 P (D) 的 相位 谱 显 示 ， 在 载波 
言 号 、 遥 测 信号 、 测 距 信号 分 布 的 主要 频段 相位 谱 有 明显 相关 条 纹 ， 其 它 频段 的 点 由 于 信和 号 
较 弱 ， 相 位 分 布 弥散 ， 条 纹 较 弱 ， 验 证 了 不 同 台 站 接收 的 探测 器 VLBI 下 行 仿真 信号 的 相关 
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图 6 仿真 两 站 信号 互 功 率 谱 图 (a) 幅度 谱 ; Cb) 相位 谱 与 条 纹 
Fig. 6 The cross-power spectrum of two station simulation signal. (a) The amplitude spectrum; 
(b) The phase spectrum and strip 


本 文 在 研究 探测 器 信号 结构 的 基础 上 , 分 析 了 包含 探测 器 下 行 点 频 信号 和 有 限 带 宽 信和 号 
的 VLBI 终端 信号 仿真 方法 ， 构 造 了 中 频 信号 仿真 VLBI 数字 终端 接收 信号 。 采 用 最 小 二 乘法 
可 构造 出 满足 精度 的 时 延 模型 ， 用 于 模拟 多 普 勒 效应 。 通 过 仿真 数字 终端 信号 与 实际 信号 的 
比较 以 及 对 不 同 台 站 的 仿真 信号 作 相 关 处 理 ， 验 证 了 探测 器 VLBI 信和 号 仿真 方法 的 可 行 性 。 
后 续 建 立 仿真 各 模块 误差 分 析 模型 ， 完 善 仿真 系统 ， 为 月 球 探测 器 VLBI 信号 仿真 的 应 用 
乃至 火星 探测 器 VLBI 信号 仿真 的 研究 打下 良好 的 基础 。 
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Simulation method of Probe VLBI Signal 


ZHU Yali!+, ZHENG Weimin!?3-", TONG Li!, TONG Fengxian!, ZHANG Juan!, LIU Lei! 


1.Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 20030; 
2.Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008; 
3.Shanghai Key Laboratory of Navigation and Positioning, Shanghai 200030; 
4.University of Chinese Academy of Sciences, Beijing 100049. 


Abstract: VLBI (Very Long Baseline Interferometry)technique is an important method for 
orbit measurement of deep space probe. In order to research the data processing methods of the 
VLBI center and evaluate the measurement capability of the VLBI system, at present the routine 
method is observing actual deep space probe by carrying out large observation experiment to get 
probe signal that received by VLBI terminals. Compared with actual observation, simulating 


probe signal based on the characteristics of the track and beacon is a better choice. The actual 
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downlink signal is radio-frequency signal, which will be converted to intermediate-frequency 
signal, then collected by receiver. In simulation, we construct the intermediate-frequency signal 
directly to reduce the amount of calculation, and generate the digital terminal signal via 
resampling. Because of the relative motion between the probe and the VLBI station, the frequency 
of the signal shifts accordingly, which is called Doppler effect. To simulate the actual probe signal, 
we build the delay model, andinvestigate the simulation methodabout theDoppler effect of the 
typical point frequency signal and finite bandwidth signal. By comparing the simulated signal with 
the actual signal received by VLBI terminal, and analyzing the VLBI correlation result, the 
feasibility of the simulation method is verified, which paves the road forthe following VLBI 
application in the moon and Mars. 


Key words:signal simulation; deep space exploration; VLBI;beacon 


